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ABSTRACT 


Diffraction radiation is that electromagnetic energy 
which is caused by a relativistic charged particle passing 
through an aperture in an opaque material. Ter-Mikaelian 
solved for the diffraction radiation from a point charge. 
This paper discusses the phenomena resulting from finite, 
relativistic charge bunches. 

Using the Huygens-Fresnel principle, diffraction 
patterns from spherical and cylindrical charge distributions 
are found and plotted. For charge bunch sizes less than the 
radiation wavelength, the results are almost identical to 
those for point charges. 

The radiation pattern is composed of two regions. The 
"transition region" is characterized by a strong peak at 
9 = y 1, the Lorentz factor. The "diffraction region" 
consists of a series of peaks and nulls in field strength 


typical of the standard plane wave diffraction pattern. 
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A. BACKGROUND 

Pitpseeteimacdalanion OCCUrS when a relativistic charged 
particle passes through a hole in a conducive or dielectric 
material. This phenomenon is very well understood for single, 
infinitesimal charges. Ter-Mikaelian computed the diffraction 
radiation for a single, relativistic electron. (eer whi eeLess 
understood are the effects of charge bunching. Presumably, 
bunch sizes on the order of the radiation wavelength will 
Significantly affect the diffraction radiation intensity by 
destructive interference. This paper discusses diffraction 
radiation from charge bunches with radii greater than zero but 
less than the radiation wavelength. Diffraction radiation 
will be computed for spherical and cylindrical bunch shapes. 
In addition, the effects of bunch size, beam energy, and 
aperture radius will be explored. 

When electrons are travelling at relativistic velocities, 
they can be treated as "pseudo-photons." A method very simi- 
lar to that used in physical optics based on the Huygens 
Pencil emismenensused tO find the diffraction radiation pat- 
tern. [Ref. l:p. 376] Panofsky and Phillips describe the 


'Viriilaleph@temeecenmcept' in Reference 2. 


B. THE ViIREGAP we HOron 

The electric field from a "fast electron” sapcroacne. 
that of a plane wave as the electron s velocity appreder 
c, the speed of light in avacuum. Representations of sen 


E-fields of electrons at various velocities are shown below: 
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For a uniformly moving electron, one can solve the 
Lienard-Wiechert potentials to find the field as viewed by 
a Stationary observer. One finds the component of the elec- 
tric field perpendicular to the electron'sS motion 1S veryenien 


greater than the parallel component and 1S given by 





il eb 
Ea ae Y 
aL eee (b* Regs) "= 
where: 
e = the electron charge; 
b = the perpendicular distance from the 
electron Sepanit; 
y = the LorenezZ sactor — ma ag ey ek. 
B = vf/c. 


CHemwecatme ley sents stvelad as a sharp pulse of electromagnetic 
energy. The magnitude of this pulse is highly dependent on 
y, which defines the relativistic energy of the electron. 
In addition, the duration of the pulse decreases as the 
electron energy rises. The electric field of an untra- 
relativistic electron will appear very nearly as a delta 
function in the time domain. A representation of this 


electromagnetic pulse is shown below: 


The Fourier transform of E, (t) gives the frequency components 


Sretnis pulse: 


_ il 
et = a _J E, (t)e Gite 
“ye -2 00 atwt 
- Y eS nares 
8nve -m — [(vt)* +y yo a 
ete © 
re 1 a 
_ e e Ee 
7 By 2 


where: 


valae 
Eyres 
Integrating, Panofsky and Phillips found 


e WD), (we 


dye = na res 
7S ane by Me al 


where Ky is the modified Bessel funetion of the second Kine 


K,.(x) approaches infinity as 1/x near x = 0 and goes to zero 


1 | 
as ete for arguments greater than one. Panofisky and Phiiseigeee 
as well as Ter-Mikaelian in his development of diffraction 


radiation, approximated Eo by equating 


(2 for o < 


Tie o 
0 for ee 


In his treatment of diffraction radiation fvom a Saaeee 
particle, Ter-Mikaelian assumes that the Fourier components 
of EL are spatially limited to a circle Of radius Bb ~ Ave 
He states that the particle will radiate at a frequency w 


only when iy > Lor where ry is the characteristic dimension 


Of the aperture. “|Ret 22377) 


Cx DIFFRACTION THEORY BASED ON HUYGENS ee RINCreae 
Ter-Mikaelian used Huygens' principle to compute the 


diffraction radiation pattern from a charged partr1elle passcmag 


he 


Piroueg@ a Cpooimmicqwmtie ENOught Of eCach point in the aperture 
as a source of secondary waves and then integrated over the 
aperture to find the total field at some distant observation 
point. The field at each point within the aperture is 
assumed to be unaffected by the obstruction. In practice, 
edge effects will be apparent for small apertures. 

Huygens Principle states that: 

every point on a primary wavefront serves as the source 
of spherical secondary wavelets such that the primary 
wavefront at some later time is in the envelope of these 
wavelets. Moreover, the wavelets advance with a speed 
and frequency equal to that of the primary wave at each 
point in space. eee pps OU o | 

Using this theory, one can describe the propagation of 
a wave past an obstruction or aperture by mathematically 
summing the contributions of each "Huygens wavelet." Fresnel 
expanded upon Huygens’ principle by stating: 

every UunOoSserMeeead point Of a Wavertront, at a given instant 
in time, serves as a source of spherical secondary wavelets 
(of the same frequency as the primary wave). The ampli- 
tude of the optical field at any point beyond is the 
Superposition of all of these wavelets (considering their 
amplitudes and relative phases). Reon O35 30)] 

There are some problems in using the Huygens-Fresnel 
principle, most notably the effects of electron oscillators 
at the edge of the aperture. Remember that Ter-Mikaelian 
assumed that the field within the aperture was unaffected by 
igewousmrnc tome OSC lating electron “clouds at the edge 
of the aperture will produce their own contributions to the 


Motte lect mice kinwithineeme aperture. The field is then 


made up of a part from the incident electric wave and one from 


itel 


the induced oscillations at the edge of thewopctrver 1 
However, for apertures larger than the wavelength of the 
incident radiation, Huygens-Bresnel Qitfiraeeren  enieom mee 


expected to provide quite accurate results. 


D.  FRESNEL AND FPRAUNHOF BRVPi Pe EAC ren 

If the observation plane is located relatively far from 
the aperture, one can ignore the vector aspects of wave 
propagation and use the scalar Huygens-Fresnel principles. 
For each small source at the primary wavefront, the contribu- 


tion to the field at observation point P(x77) is eemveree 
X,Y) = E, (x,y) exp 1 (Wisk aids i 
The exponential term provides the phase information and 


the r term in the denominator takes into account spherical 


Spreading. A diagram of the problem is shown below: 


ee 
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In the typical diffraction problem encountered in the 
study of physical optics, a plane wave is assumed at the 
aperture, so the electric field across the aperture has no 
Spatial variation in amplitude, polarization or phase. When 
dealing with the radiation from charged particles, however, 
One must consider that the field is radial with amplitude 
changing with distance from the charge. In our problem, we 
locate the observation point on the X-axis and deal only 
with the x-component of the field at the aperture. The field 
is symmetric about the Z-axis. 

The total field at P(X,0) is given by the superposition 
of all the Huygens wavelets from all the sources at P(x,y) 


of size ds across the aperture. 


ebqo sl, pea) AG 


Be (X,0) = J E (x,y) - 


Fx AOC EULE 


iiss a yery Gifficult integral to evaluate in that the 
exponential term has a modulus of 27, making it difficult to 
make an asymptotic approximation. 

The phase factor can be approximated by taking the binomial 
Sapansltoneeor on and approximating it to the first order. 


2 
oe = ae + exc 47 =) = hee 0 ie = 74).0 aes x? + y 


R-(1 - 2Xx/R? + see + pie) 


i 


where: 


Bee ez 
SO 

r= RL © 2Xx/R° | (2 
The binomial expansion 1S given by 

GU et ax + —— ae 
TO f1iESst Order, 

Y = Ree hee ee ee 


The total field at P(X,0) becomes 


qi (wt-KR) { E ax (XY) 


ApereEuine i 


Hex 1/2 SRR) 


ae ds 


For large R, one can replace the r in the denominator with 


R and the integral becomes 





qi (wt-KR) -—ik ( 
Pr Xx Iss YX 
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If Ris much greater than the aperture size and the 
phase difference from all of the incremental sources, ds, 
across the aperture is negligible, then the k (x“+y“) /2R 
term in the exponential is very small and may be neglected. 
In general, if this term is less than one radian, then the 
integral can be approximated by 

ot (wt-KR) 
a eee = ——g J Egy (ry) exp (LkXx/R) ds 
In the Fraunhofer approximation the electric field ata 


distant point is just the two-dimensional Fourier transform 


of the electric field in the aperture. 


dS) 


II. COMPUTATION OF PIE PRACTION RADI AIwe 


A. SOLUTION FOR A POINT GHARGE 
A single relativistic electron has a Fourier spectrum, 


as derived earlier, given by 


2 eee 
1s = [ cs )K 


b 
- ae 
An EGPVv 


1 ‘vy 


Using the Fraunhofer method, the field at a distant observa- 


E1On POINT is Givens, 


ie 
_ e se 
E x (%19,2) rom a f Ea,x'xryle 


7 ads 
Aperture 


where the x-component of the field at the aperture is given 


by 


E = —y2 1 (*2)K (2) cos 6 


Changing to polar coordinates gives us 





X 
1 (wt-«R) O 271 
ae a C2) fore 9) ep ato tts 
0 0 41% by Mi Soe 
“a 
1 (wt-«R) O 27 
ae ; = f Py x, { aoe 9 @t (KbX/R) cos 6 dé db 
Ant EQVv O be vy 
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Integrating over dé yields 


: A 
1 (wt-KR) O 
e earn wh kX 
ei (wt-«R) 1e Ww % x x me 
7 7 Cie wy %2, (a2 wy “tem! Sil 
VY 
Ww x Ee 
~ 55 Jy (kA, RIK (2) 


miv~ch contains the final form of the diffraction radiation 
pattern. As Ter-Mikaelian did in Reference 1, the radiation 
aioe FOUN woyectiomracting the particle field from the 
results above. The particle field is found by letting Ne 
go to infinity in the integral over the aperture. It is 
interesting to note that as A goes to iced eee the ed it fiece— 
tion radiation field strength goes to zero. Physically, at 
this point the electron no longer "feels" the presence of the 
obstruction. 

Ter-Mikaelian, in Reference 1, computes the diffraction 
radiation from a charged particle passing through a circular 
aperture. Using the method outlined in the previous section, 


Be dectrermines (tor OAL ell) 


IE, 


where: 


UJ 
a = =; 
ra 
q = kK sin 6 = the projection of the wave veceem 
k on the plane 2 ==0; 
1s the aperture radius. 


This result appears to be significantly different from the 
expression derived earlier. In reality, however, the two 


plots are virtually identreat (See Appendix B.) 


B. SOLUTION FOR A SPHERICAL CHARGE BUNCH 

The diffraction problem has been solved for single 
charged particle wadracion, 

The process for a finite charge bunch can be solved, as 
a first approximation, by assuming a spherical charge dis- 
CHhIbuULIOnN OF Rradaus a, with constant charge density, 0. 

In the frame of the moving charge, the electric field 


is radial with a magnitude given by 


< 

z2 EG for wees, an 
O 
E = 3 
ne 

Oey ees ay 

3e€ 24 
This result is easily obtained using Gauss" law. A repre- 


sentation of E(r) 1S showne below. 
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meer rnd the field from the moving Charge distribution in 
Picertrame Of a Stationary observer, one must find the Lorenz 
EeaiislLormation Of the radial field. With a spherical charge 
distribution, an observer views a field emanating from the 
eemecr Of the Charge bunch; the field from the uniform 
spherical distribution appears to be a point charge equal to 
the total charge of the sphere. If the observer is located 
within the sphere, he still "sees" a point source, but of a 
reduced strength (that portion of the sphere's charge enclosed 
by a Gaussian surface). Conseauently, the Lorenz transfor- 
fattomeot the wmoving Spherical charge bunch must nave a form 
Seebar tO that ©f the point charge. In the lab frame, one 


finds a perpendicular electric field of 


4/3 Tap ~2) 
Tie eee acts) 2 ey 
2 Oo [(vt)”~ +y “b’] 
ott) = 
| 
3 -2 
av eis Cp 
a — Omens oa 
1 (vt)? ay 2471372 O 


ile 


This form Of E (t) 18 Only aeppmee<imatecwnew som a,j: For 
small bunches (a. <e me ar the error introduced by this 
approximation is negligible. 

Since E (t) of the sphere has the same time dependence 
as that of the point charge, it must have Fourier components 


of the same form. 


4/37a_o 
522) x, (2) for b> a, 
4n“e bv ! Y 
oe : 
afore 8 (2) KN] coi ae 
4t EQPv Y 


Again, Ek is negligible for X < b/y, we make the additional 
requirement that the wavelength be large compared to the 
bunch size in order that the bunch radiate coherently. 

Again, one proceeds uSing the Huygens-Fresnel principle 
to find the radiation pattern from the spherical bunch passing 
through a circular aperture. With the observation point 
located on the X-axis, one expects to see only an X-component 
to the perpendicular field. Due to symmetry, only the x- 
component of the field across the aperture contributes to 
the total field at P(X,0). The transverse electric field 


at P(X,0) is given by 


ot (wt-KR) 

E Cx 0 2,) SO f lo (xe 
roe Zo a 
(POV SHOE BUCS 
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ieee tee x-cOlMpomene Of Che field at the aperture is given 


by 
6: 
4/37a_o 
5 = [(22)K, W°) Icos Cmeror bb < ao 
4n-e bv Y Y 
R = 
a,X 
' 3 
oi one 21 (G2) K, #2) Icos COte ie aa 
at EQRv Y if 


Changing to polar coordinates yields 


1(wt-KR) ; 
_ e 1(kXbcos 6/R) 
EDX 7 ——z— JE, y(b,8)e ds 
where: 
as — oe arcls db 
Finally, one obtains 
1 (wt-«R) 
e i (kxbcos 6/R) 
eee b dédb 
1 (wt-kR) we Pa Ss (kx 
= ee i f asa 8 M2) K, (2) Icos get! S/R), ae db 
0 0 4n%e by er. | 
‘o au 4/3na5p wh wh 1(kxXbcos 6/R) 
fe 5 [ () K, (=) Ieos Be b dédb 
a, 0 4m € bv 0 ad 


Zagh 


We shall first solve for the integral roen >= 0m oe 


° 
f 


we have 
“5 ei 3 
f f a (C2)K, (2) Icos Ge7 kxbcos 6/R) b dédb 
0 0 4m ebv at Vi 
4 3 or . 
eaves f b eK {cos ‘ca Sees cael. 
4m -e_v 0 vy ai 


O 


The integral over d6@ is equal to (271) J, (KXb/R) . The integral 


becomes 


As an example, one can expect to see parameters in an 
experiment on diffraction radiation in the micro-wave regia 


as follows: 


f = 8 GHz 
a = 0.5 cm 
Oo 
Ymin > AM 
Vile Moe ae” cm/s 
-l 
k= 217k?) Sale 5c 
z c > 
X/R = Sin 6 where 6 is the angle between k and the 
Z-axis. 


Ze 


Poiigmenecsesvalulcs »mene Maximum values of the arguments in 


K, (x) and Jy (x) Mimeie=IiLeqralerrom 0 £o a, are 


Shy le. LTC 


(wb/yv) 


1 


(kXb/R) 8.38 x10 


With these values for arguments, the two Bessel functions in 


the integral can be very closely approximated by 


K, (x) = IL ss 
and 


Bo) 2 


(x) 


The integral becomes 





a 5 
_ 4/3n9 ~ 2° .31,kxb | DR ais ae 
Zi Ane vy { le) x (a) db = Za =e kK) (=) Simo 
0 4m EV 
O 
We next evaluate the integral from a, to 4): Since the 


integral here has the same §-dependence as in the previous 


case, the integral becomes 
3 n 
4/317a_p O 
2 (wb wh KXD) db 


271 —~— i ) le 
Ante rs ee aves | RR 
O O 
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Note that in this case, we have moe made Ene aporoxiiaeeems 
for the two Bessel functions since b will in practice extend 
to values much larger than one wavelength. This integral 


can be easily solved and yields 





O 
wb,., ,wbo,. ,kXb 7 1 
(PAGS = Cae 
40 R VY 

KX Wt KX i KX, Wt 

«gh Cy Ca) Re a oe eae 
a (OK, CaP), ed) = Bg (eK, (NH 

OR Lb wy a aes SnD 3< D VY 


Combining the solutions of the two integrals and sub- 
tracting the particle field, the total transverse field at 


P(X,0) is given by 


3 Z 
1 (wt-KR) 4/3Ta_o a 
a O i. kX ee W 1 
Oy eet 0) Z (2 Fama Glee, (=) wR? , aye 
aie YV 
kx ears O W noida Xs 
l1 ot (R) Ky ye 2 ' R i Yv J, R )Ky Ca 
-a {(S)k Pua FO | mt Oy (2) 37} 
OR lL” ¥V9s2 ae vo Re ae 
(esenay? «(Be 
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e. PPP TOS rABLAITION FROM A LARGE SPHERE OF CHARGE 
POu adtrfraction problems where the radiation contribution 
Monge ene chamge bunch 1S Sligmificant (i.e., ae & Ons 


one can not make the approximation used in the previous section 


where E, had the same time dependence for b < a, as iki@aralael foes 


i 
b > a,: Mimegicuseemion, we shadl gain a more exact under- 
Peano ngeorecenesproblem by determining the diffraction radia- 


iron £rom a large sphere of charge. 


Using Gauss's law, one finds 


3 
eae sdhcibe b BOr JG fea 
4une a —= "© 
ne 
B ss « 
A 3 
0, 47 302 b a 
Ane 7 ia — Oo 


where b/r provides the perpendicular component of E, as shown 


below: 
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This relatilon@recuces 26 


Aj 3ana-b 
Oo 
3 LOR Wee! -a 
Oo 
p : 
ee ~ Ate 
4/3 1b fOm = a 


fh GeprEeCSen eat le itor EY is shown below for several values of 


b: 


Small Sphere 


App vex mation tor bda, Large Sphere 


Approximation for bt a. 





Figure 1 


TNS LUNCELeh stor som a, can be integrated numerical 
only. However, one can make a very close approximation for 


E with b < ay bY wWeat Bog 


26 


4/31mb(1.5a_-Z) 
ae es FO paca a7, < T.5a. 


O 
L.5a_ - a’—b* 
Fo = Ny 4/31b 1 Ose = Ly, oa 
L Ate O O 
4/3mb (1. 5a_+Z) 
————-. for -1.5a <2 < - Vat -b* 
Stl 5a - Va2-p? 
O O 
This approximation is shown on Figure l. For most values 


of b, this linear aoproximation to the Lie dependence of 
EY (r) outside the sphere lies very close to the actual curve. 
This is particularly ttrue as b approaches ays where most of 


the diffraction radiation from within the sphere is expected 


meamorlginate. Replacing zZ with -vt, E (t) is equal to 
4/3naep 
n@ lage. > acl 
(b+ (vt) 7) °/7 
ay stip ( 5a _+vt) 
——E———————— for b<a DY ac=b- Gi= Vier <i eal 
O 1.5a_ - WS -b 
eS O O 


= 
4/31b jig « Jo) =< aor 7 Vat—b* <= < Vat -b" 


AV SU on ANE Sie) aoNany 
Se noiemuoe mci = log. < —VE < = ce je 
eee) O O O 


1.5a_- Wes 
O O 


Rema olaetyiotac electron bunch, the Lorenz transforma- 


tion gives 
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4/3Ta_b 
5 a 23/2 LOR mon a 
(b°+(yvt)“) . 
4/3nb (1. 5a +yvt) Vat —p* 15a 


for .b <a 


eee O Y Y 
1.5a, - Vat-b* 























E(t) = ne ss 
Ate ~ V92-p2 Va2-p? 
4/3tb fOr; %s < az = < -vt < = 
O x ¥ 
4/37b (1.5a -yvt) ~1.5a, -Var—p 
a Tom) 0S = ou < -vt < 
DS O a Y 
1.5a - a_-b 
Oo O 
Or 
4/3na~b 
eae rere for bra 
2 ? 
(b+(yvt) ye = 
4/3mb(1.5a +yvt) ~1.5a, - Var? 
Boje Joy << aor vy <ies- wy 
Bit) = oa Vae—b* 
Ate Vae—p" 
4/3nb for b<a, |t| < ———— 
Oo Vi 
4/3mb (1.5a —yvt) Va‘ —b* 1.5a, 
——— for b < a_, ———- < t < ——- .. 
; 75 Oo Wey VY 
5a la a -b 
O Oo 


The diffraction radiation found by integrating over the 
aperture from a, to *%, for the large sphere of charge will 
be identical to that found in the small sphere approximation 


and need not be discussed here. 
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Ue msceutenmwd ileal scisseonly that radiation originating 
from the aperture for 0 < b < a,: For Chis range of b, the 
aiiaetneonpotentseor the field are found by writing E, (t) 
as a difference of two simple functions whose Fourier trans- 
memms are tabulated in any introductory text on Fourier 
Mmiimiyvsisc. EF, (tj) 12S written as the difference of two tri- 


al 


angular pulses, as shown below: 





It is easily seen that 





1.5a, 
KE 7 \, etl a) Sy way? 
wh 4ne€ 2 
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4/3mb (1.5a_) Sin’ ( 2vy ) 
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A diagram of @ typical contribu ronmtomn from each function 


1s shown below: 


ar 
2 / ee 
a SIN \4 9 ae 
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One can easily see that Pet is cut off at high frequencies. 
The field at a distant observation site P Oe.0 . 7) age 


given in the usual manner by 


i(wt-eR) %o 9 2a 
BE, y(X0,2) = = —— ff Ecos get KXbcos O/R),. a, 
Dr X Z 0 0 Wd 
: a 
1 (wt—KR) O 
= kXb 
: 2ni Wd, aa) ole 


This integral can be evaluated numerically to find the dames 
fraction radiation pattern from a large Sphere of (chamger 
When this is done, the results are not accurate enough to 


make detailed predictions regarding the diffraction patrons 
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They do show, however, that the diffraction pattern from a 
large sphere of charge is consistent with that found for 
icmmounehes. ine field strengths are significantly less, 


but the pattern of peaks and nulls is roughly the same. 


D. DIFFRACTION RADIATION FROM A CYLINDRICAL CHARGE BUNCH 
Peete G remeste ma eotatloOnadmyvecharge distribution with 


axes as shown below: 





{= buach len th 
Qs bunch radius 


We wish to find the electric field at an observation point 
Pec by cylindrical symmetry, the electric field at an 
observation point on the X-axis will have only : and k com- 
ponents. We are interested primarily in the i component 

EInece, Peete thew relativistic electron» bunch, Ea Ey 
Eneorader eo solve for EB}, , one must fipst tind the electro- 


Settle soomutlOmmror a Cylindrical bunch. One can make an 


Hop EO iat lO spy meatincg tne cylindricalsbunch as a line 


Sil 


charge of linear charge dense (1a 





s”6~—Cr FO <a ee 
O O 
linear charge density 1S given by oe A diagram of the 
problem is shown below: 
x 
P(X, 2) 
r- 2’ 
2” AE 
FOG =D > wa 
oTa . 
DO A ZA), = GZ 
4te€ (x? + (Z z1)*) 74 
Let u = Z-z; du = dz 
nal 7 Ta Z 
Ea = "oe 2 du vio a | J 
L 4te o-L (p24) 2/2 4te 2 (p74 17) 1/2 -L 
mac 
ao Z Z-L, 
4ne b (b742.7) 1/2 (b7+(2-L)2)2/2 
For a slow moving bunch, located at the aperture at — = 9, 
one can substitute (-vt) 


for Z and write the time varyiug¢ 
field at the aperture as 


B, (t) 7 QO, (-vt) vt+L 


___ (Vb) a 
(bo eGo eo ete 


oe 


Poe lativistleschnarge bunch, the component of the 


electric field perpendicular to the bunch path is given by 





pmb (vty vty) 
e272 Fy) 2) 172 
Ae). = 
oma 
ae a 7 * oh ee UP) 
4meb (b+ (vty) ~) / (b°+(vty+L) ~) / 


Ree. < a, 


NRO) Gag OG Yaa) 
O 


As in the previous section, the Fourier components of 


im t£ield for b > a, are given by 


el > Iwt 
ee >. | OC (te at 


wt 00 lwt 


an Ane Oe Wea vty4L) 2) 2 





BEyelUating the first integral, one finds 


00 1wt eo vtyt+L 
{ (vty+L)e te = | a, 


Vb 4 (vty4b) “ my [,, (CEYHLy 
' b 


it 
e 
dt 


Peeing ey tiniyb, tt — (bi/vy)—-(L/vy); and dp = 


the integral becomes 
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(avy On) Ges 


; b 5; Wb 
© iw (— y-—) WL i—yV 
b VY " vy — co VY 
vy S = ay = a eX f Sa 
1 eae 
Using the equation Ql? = cos 6 + isin §and remembering the 


properties of even and odd Fourier integrals, the integral 


1s rewritten as 


Using the same techniques to evaluate the second integral 


in £ ..,. One finds 
Wh 








2n 4neb 


omar oe! 
a Ome 42) wb VY _ 
‘s amen * WS Sp? ‘ I) for b> a, 
E = 
Wh 
i onb? 2b. ub, Sy 
7 Sr Fy" e = AL) ore 16) a, 


As previously explained, the Hankel transform of the field 
at the circular aperture is evaluated to find the field at 


observation point P(X,0). 
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a WL 
1(wt-KR) . O 27 Z 1— 
= 2 } kXbcosé 
E = | f f ollie eon EAE VY-1)cos6e" = / ape 
0 


Simplifying, the equation becomes 


ei (wt-KR) 


al 
Cy ip 
x] f grater oe 4 cos dbs 
te wey 
O 0 


O Zi KX] 
+ at f if bK (2) 2 = EOE 8 db a 
Ca 0 a 


As in the problem for a spherical charge bunch, integrating 


Over d@ gives 


1 (wt-KR) a: ye 
E(x,0) = =—~—(- Tas? NETS 
a a 
3., ,Whb kXb 2 wh kXb 
O 


3 


Evaluating the integral from 0 to aj, one makes the usual 


approximation 
Fo 
3 wh kXb ~ VY 2 kXb 
f b°K, (7) 5 (=e) db = Seno) eee 
Since in this range, K, (x) 1/x. One finds 
ay 5 KXa , 
YY ff ps, (2 yap. = Se ee 
W 9 i ees W 


Next, evaluate the integral from a, to toe 











AL 
2 2 wb kXb 2 1 
3) ; if OE ae eee ao EX 2 ey? 
O R VY 
Kove WAL KXA uy7t 
kX O O W O O 
x (A tJ 5 | =) LS are = vy 1! R )K, err © 
KkXa Wa kXa Wa 
kX O oO W O O 
a tea. R PS a ay | R Ro (a) 


Since X/R = sin 8, the total field at the observation plane 


ale 
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The effect of bunch length is in the (e vY-1) term. ox Oelig (ol= 


ing the exponential in this term gives an imaginary and 


real part of the solution. 


ecole WL 
a (—) f2 
(=) =e 2i a7 
WL 
= Wale vy /* stim (Ss 
2VY 


Si, 


At first glance, this function seems to give the diffraction 
an oscillating dependence on radiation frequency for all 
frequencies. One must remember, however, that the Kj (wa /vy) 
term in the solution very nearly cuts off the radiation at 
> ae The real and imaginary terms in this expression 
give an additional phase factor to the overall solution for 
diffraction radiation. For wavelengths much larger than 
bunch size, the real term is approximately zero and the 


dependence on bunch length is given by reo 


Mi 


E. APPROXIMATIONS AND LIMITATIONS 
The reader iS cautioned to remember that in making the 


Fraunhofer approximation to the diffraction problem, one 


assumes: 
1) The diffraction angle is small (< 1 rad). 
2) The phase difference from all "sources," ds, across 


the aperture is small (k (x? +y2) /2R) ares 


3) The aperture size is greater than the wavelength of 
the incident radiation. (Neglect edge effects.) 

4) Since charge was assumed to be uniformly distributed 
throughout the bunch, the quantized nature of charge 
was ignored. This can be done if the electron den- 


Sity is sufficiently high (n~173 < }). 


5) Charge bunch sizes are assumed to be less than the 
radiation wavelength to ensure coherent radiation of 
the bunch. Larger bunches have smaller intensities 
due to destructive interference. 
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teers AFFECTING DIFFRACTION RADIATION 


Sees Gee 


imei lewLelcmchapter,  Girtrackion radiation was 
computed and, for a given frequency, was found to be depen- 
dent upon several parameters, namely y, bunch radius, aor 
ema aperture radius, Noe In addition to these effects, the 
madiation IS Signiticantly cut-off at high frequencies. 

Remembering that both the small sphere and cylindrical 
charge bunches had Fourier components proportional to 
(wo/vy) K(wb/vy), one can make the approximation that the 
melds are small for w > vy/b. (When the argument of Ky (x) 
Pome (xX) = .F0l9 This gives a field intensity of about 
One third the maximum value.) If the maximum value of b is 
the bunch radius, ays then the maximum value of w is given 
by w = vy/a. > One expects radiation at 20 ax to be 10 dB 
less than the peak values. 

Figure 2 can be used as a rough guide to the frequency 


dependence of observed radiation. 


X 


Figure 2 


Se, 


For the large sphere, the Fourier components of the 


fields were found to drop off as 


One can say that, the filelaers cue ore aioe vy/a.- 

It is easily seen that, for each bunch shape discussed, 
the diffraction radiation is very near the peak values at 
microwave frequencies. A previous experiment at the NPS 
Linac discovered diffraction radiation in the microwave 
region (f = 8 GHz) from a circular aperture [Ref. 4]. For 
f = 8 GHz, wbh/vyK, (wb/vy) 1s very nearly equal to unity. 

The plots of diffraction radiation vs. 9, the off-axis 
angle, consist of two distinct regions. The first is charac= 
terized by a strong peak at 6 = ae Since the peak at 
§ = va 1s typical of transition radiation, this regione 
the curve is called the "transition region." 

he ae eo the radiation pattern appears more as the 


typical diffraction pattern of peaks and nulls. This region 


is called the “dititraceien weqren. 


A. DEPENDENCE OF DIFFRACTION RADIATION ON y 

The y-dependent portion of the function describing diti- 
fraction radiation tends to be very complicated and hard to 
understand analytically. The strongest dependence is in the 
transition region, where a very strong peak in field strength 


occurs at 6 = a? In this region, the field strength is 
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wervyeneariy Given by (k sind)/((k Sa HAaaea Ne. The 
ZaeUjeimehmr Nts upeak TS directly proportional to y. In 
the diffraction region, the effect is less pronounced. 


eMewexPpects £O See more diffraction radiation energy from 


electrons of higher energy. Figure 4 shows the dependence 
Werner ackton radiation on y. The curves for a cylindrical 
Pumciewitl be identical to this figure. Increasing y tends 


Eemocharpem the diffraction pattern, increasing the strength 
of the peaks in field strength. 

For pseudo-photons of a given frequency, the transverse 
melds tena to extend farther out for increasing y. Alterna- 
tively, the field for a given frequency and radial distance, 
b, tends to be higher for an increasing y. Figure 3 demon- 
strates this point. If one compares the same Fourier com- 
ponents of two relativistic electrons with different energies, 
the more energetic electron's fields extend much further out 


than those of the less energetic electron. 


(aby K (we) 





bacuae 3 
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ea tiie dmEaeremmined that, for a point parere le 


with Wt SVY << 1, the y-dependence was given by 


osu ieee 


; 2 ieee 
(k sin®é) + 


Rie, -cependence found in the last Chapter for both cylin- 
erececanicespaerreal charge bunches 1s consistent with that 
eeler—-Mikaelian. One expects the fields to go to zero for 


Meie-relativistic particles. 


B. DEPENDENCE OF DIFFRACTION RADIATION ON BUNCH RADIUS 
Throughout this paper, we have assumed that the bunch 
radius is smaller than the wavelength of the radiation. When 
this is the case, one may assume coherent radiation from 

the bunch. 

By examining the equations for diffraction radiation for 
mimemenree Dunch types, one can see that diffraction radiation 
tends to decrease for bunches with larger radii. For small 
bunch radii, a general decrease in the field strengths 
from the peak levels for infinitesimal bunch sizes is observed, 
as shown in Figure 5. For cylindrical bunches of equal 
Seager wene Gdititaction radiation field strength is inversely 
proportional to L, where L is the bunch length. This is 
expected since the radial fields at the aperture are propor- 
tional to the charge/length within the bunch. The ODN 
factor also specifies a length at which coherent radiation 


from the bunch can be expected. 
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One can easily understand the effect of bunch size by 
examining the static radial fields from a sphere of Glictece . 


it omecterrom ascylindrical Gharge bunch will be similar. 


= 


Figure 6 


Figure 6 compares the fields from two spheres of radii 
a and as: As shown earlier, for an infinitesimal sphere, 
Emie Llelds increase assymptotically to infinity as r goes 
to zero. For a sphere of radius aor the field increases 
with the same 1/r dependence as r decreases from infinity 
tO ao and then drops off linearly to 0 at r = 0. One can 
easily see from the figure that the sphere of smaller radius 
is the source of more electrostatic energy than the larger 


Sphere. As the sphere gets smaller, the total electrostatic 


Prmemgyeraptdly approaches that Of a point charge. 


Ce lLeEPENDENGE Or DIFFRACTION RADIATION ON APERTURE RADIUS 
A plane wave incident upon a circular aperture results 
Peete riety Aitiraction pattern... The large central 


HoxdMuUMelieintensity 1S known as the Airy disk. The 
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diffraction radiation pattern from a relativisetc cleocotron 
beam is very similar, except Enat the Ccenival ma aigomeae 
replaced by a null. This central null results from the 
radial nature of the electron fields. In the Fraunhofer 
diffraction method, integrating over the aperture for dé 
results ina Oth order Bessel function of the first kind 
for the plane wave and a lst order Bessel function for the 
electron diffraction radiation, thus explaining the Prine. 
pal difference in the two radiation patterns. 

Figure 7 shows the diffraction radiation pattern for 
cylindrical and small spherical charge bunches passing 
through circular apertures of various radii. Changing the 


aperture size affects the diffraction radiation patlecuneen 


both the transition and diffraction regions. One can define 
° Fa ZN 

a dimensionless parameter, KAOe where k = =k = the wave 

vector and _ is the aperture radius. This parameter alone 


determines the location of the peaks in the diffraction region. 
When KA, increases, the peaks gain in strength and move in 
closer to the 6 = 0 axis. Decreasing KAS spreads the dif- 
fraction pattern out over wider angles. This is predicted 
Since the diffraction radiation pattern is caused principally 


by the J (ki sin 6) terms in the solwveien. ) Note ene 


Depa 
patterns with a large number of maxima and minima are found 
only for aperture sizes much larger than the wavelength of 


the observed radiation. This is to be expected since, for 


the normal diffraction of a plane wave, a small aperture acts 
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almost aS an isotropic point source Of 4admation. “Siamede 
case of diffraction radiation, we still see the null at 

86 = 0 due to the radial nature of the fields. Since the 
dimensionless parameter is given by KAO one Can change tice 
appearance of the diffraction pattern by changing i, the 
physical size of the aperture, of k, which determines the 
wavelength of the radiation. 

In the transition region, changing the parameter, KA Os 
affects the strength of the peak at 6 = ra As ce goes 
to zero the peak at sons increases, asymptotically approaehum, 
the value characteristic of transition radiation. As KA, 
increases, the strength of this peak decreases, until, as 
KA, goes to infinity, the peak goes to zero. At this point, 
the charge no longer senses the obstruction and no energy is 
radiated. 

One should also notice that the diffraction radiation 
pattern for spherical and cylindrical charge bunches are 
identical for bunch lengths less than the radiation wave- 
length. For very large charge bunches, interference effects 
are expected to reduce the overall intensity of radiation/ 
(bunch charge), but the location of maxima and minima will 
not change. The locations of the radiation pattern's maxima 
and minima are principally a function of radiation wave- 
length and aperture size, as it is in the normal plane wave 


GLEeraction PEOb len. 


48 


TV. DIFFRACTION RADIATION POWER 


It has been suggested that diffraction radiation may 
me used for electron beam targetting or in a non-interceptive 
beam diagnostic device. [Ref. 5] 

Present diagnostic systems are based on the beam passing 
BiEaugm a thin folk and using transition or Cerenkov radia- 
tion to analyze beam characteristics. The problem with this 
scheme is that the beam is disturbed by the foil. If the 
beam passes through an aperture, one can, in theory, obtain 
the same diagnostic information without disturbing the beam. 
Bemetrrraction radyation is to be used as a diagnostic or 
targetting means, one must gain an appreciation for the power 
levels expected from a typical device. 


For an electromagnetic wave given by 


3 ei (wt-KR) 
Oo 


Ch 
II 


f 


the power per unit area, or irradiance, is given by 


Gnewcan compute the diffraction radiation irradiance by 
using this expression together with the equations derived 


mm Chapter i1. 
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As an example, the diffraction radiation power 
delivered by a 1000 amp, 50 MeV electron beam passing 
through a 20 cm diameter circular aperture will be maximized 
at 6 = 10. If one assumes 10 cm long cylindrical bunches 
of radius .1 cm and evaluates the expression for diffraction 


radiation from a cylindrical charge bunch, the irradiance 


at Z = 1m from the aperture 1s egual to 
CE 

_ Oa 2 

TI = or ES = 500 KW/m 


The power per unit solid angle is 500 KW/SR. 
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Vee CONGCHUS TONS 


PMimeiicmpaper yr tne CCuUations regarding diffraction 
Boclbalion £rom a relativistic charge distribution passing 
through an aperture have been derived from the Huygens- 
Fraunhofer treatment of the fields across the aperture. The 
calculations are very similar to the familiar problem of 
plane wave diffraction solved in every physical optics 
class, with the complications that the fields are radial 
with varying field strength in space. 

While a plane wave incident on a circular aperture causes 
a diffraction pattern with a large, central peak known as 
meme Airy disk, the diffraction radiation pattern produced 
by electrons is characterized by a null at 6 = O and is 
highly dependent upon aperture size. In addition, a strong 
peak at 6 = ve is present, characteristic of transition 
maciadteion. If Ehe aperture size 1S very small, the structure 
in the diffraction region disappears and the radiation 
pattern approaches that of pure tranSition radiation. The 
only radiation present is the strong peak at 6 = eo. Iie 


creasing the aperture size, or more properly, the dimensionless 


parameter KAO decreases the peak at 6 = vie Bice ecOmpresses 
BiewSEnUGEUhe dma Ene GQiffraction region. Eventually, as Kr 
BepEeCtehes wintinity, Ene raddation goes to zero. The diffrac- 


Mion nadlatromepattern in che diffraction region exhibits 


Sk 


some dependence on electron energy but not nearly as strong 
as in the transition region. The radiation intensity in 
both the transition and diffraction regions has a minor 
dependence on bunch size. 

Ter-Mikaelian describes diffraction radiation as having 
a strong peak at 6 = / ai followed by a broad series of peaks 
and nulls dependent on J (k sin6%,). The diffraction pat- 
terns developed in this paper are primarily dependent on 


Jy (k sin6%,) and J, (k Sin 64,), but when plotted are found 


Z 
to be in complete agreement with Ter-Mikaelian's results. 
Both results are in agreement on several key aspects of 


diffraction radiation. 


1) There is a strong dependence of diffraction radiation 
On aperture size. 


2) The diffraction radiation energy goes to zero as 
aperture size goes to infinity. 


3) The diffraction radiation energy is dependent on the 
energy of the beam, described by jy. 


4) The CE ae radiation exhibits a strong peak at 
G6 = Y~ 
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EVALUATION OF LINE CHARGE APPROXIMATION 


ee et 


BORTCY IL INDRTer tl BUNCH 


As described in Chapter II, the static fields for a 


Syeindrical charge distribution were approximated by a 


modified line charge. The linear charge distribution was 
given by 

ia — ona. C7 Cite 
where: 


eyiemparical bunehe@ radius in cm 


a5) 
ll 


Gwargesdensityv ef cylindrical charge 
puncimin C/Cm 


© 
II 


When b, the off-axis distance, is less than a, the linear 


charge density was given by A = oemTb 





For the geometry shown above, the field for the line 


charge is given by 


Se) 


ee 2-b 
2/2 We ee eee 








As a check of this approximation, the fields predicted 
by the linear approximation were compared to the fields 
resulting from the actual cylinder of charge, which were 


computed by numerically integrating 


{ O Pari 1 Sees cUs-2F T37z x drde az 
((X-r cos6)°+(r sin8)”~ + (Z-Z')~) 


p 
ATE 9 6=0 Z=0 

Figure A.1.1 shows the comparison for three "long" 
cylinders. (ie as). where the observation point is outside 
of the cylinder (b > as). The linear approximation is 
expected to work best in this situation. 

As can be seen from the figure, the approximation is 
very close for the cylinders of lengths 1 cm and 2 cm and is 
identical for the longer 10 cm bunch. The radiate rerae 
become very nearly equal to that of an infinitely long line 
of charge very close to the ends of the bunch. 

The linear approximation is expected to perform worst 
for a relatively short, thick cylinder of charge when the 
observation point is on the interior of the cylinder. 

Figure Al.2 shows the comparison for these cylinders 
ranging in relative length from L = a Or nae an: For the 


shorter length, the approximate fields are found to rise 


slightly more sharply and to higher levels than the actual 
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fields. Even in this most extreme example, however, the 
“@eroximarion should be adequate for our purposes. As the 
bunch length increases, the actual and approximate fields 
fee LOunad tO agree very closely with a slightly steeper 
rise in the linear approximation. 

imeconclusion, we find that treating the cylindrical 
charge distribution as a modified line charge offers a very 
Satisfactory approximation to the actual fields present. 
DenmgeehisS approximation for a static charge bunch, one can 
Compute very closely the Fourier components of a relativistic 


charge bunch. 
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APPENDIX B 


DIFFERENCES WITH THE @HEORY OF VER-VIRAnE ea 


In Chapter III, it was brought to the reader's attention 
that the equations derived in this paper were slightly at 
odds with the results presented by Ter-Mikaelian in Reference 
1. In his paper, Ter-Mikaelian stated his results for an 
electron passing through the center of a circular aperture 


Of Ssadiis a Pou wa ./VyY << 1 {in Gaussian@units): 





ee q 

E, = ; 5 5 J, 'q1,) cos w 
ee Cl Sse! 
1e ; 

E = J (qr_)sin w 

Y ae S| toe . 2 

where: 
a = w/vy 
q = ksin6 = projection of the wave vector on 
the z = 0 plane 
y = the angle between Gq and the X-axis 


These equations give a radial field equal to 





= 1e G 
ey, = 5 5 J 


qi) 
Zia ac + : 


3 | 


Since this equation appears to be inconsistent with the 


equations derived in this paper, this appendix will explain 
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Helenoct uae method used tO Obtain the Bessel function in 
the equation. 
Hiemiinteecralwovwemeag in the Fraunhofer diffraction problem 


was written as 


ZA 
i cos 6 e 
@ 


i (kKbX/R) COS 8 44 ay 


The cos § term in this integral is due to the radial nature 


of the electric field. Using the expression cos 6 = lave 4 


2, the integral can be re-written as 


PAN 10 Sole 
{ (e = yet (kbx/R) cos 8 a6 
@) 
2m fico eeose ) (oa co sch 
ik R R 
= — e +e dé 
=) 
@) 
ie 2m 1 (mvt+ucosv) 
Since J(u) = —=— f e Aver Ox 
m 27 0 
m 2m i (mv+tucosv) dv 
oc) ee nc : 


the integral over d@ is equivalent to 
Jal _ 7 
Tid, (kXb/R) + Ti JN S873 = dy a J_y (7) 


5) 


Since J! (x) = 1a. _, (x) - J, jaya) wii = C7 eis 


expression reduces to 


kXb 
= i ( 
2 J 5b? 
which equals 
; kXb 
21a Jy (-27) 


since Jo ) - ~J, (x). 
As shown earlier, the sin 98 dependence of the diffraction 


pattern 1S given by 


il 


2 So ig ah; Ce 2 ee - ad, (qo) Ky (aro) 5] 


This appears to be quite different than that found by Ter- 


Mikaellan. 


However, if one plots the two solutions, he will find 


they are identical. 
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